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Analysis of spacecraft vehicle responses to atmospheric wind gusts during flight is 
important in the establishment of vehicle design structural requirements and operational 
capability. Typically, wind gust models can be either a spectral type determined by a 
random process having a wide range of wavelengths, or a discrete type having a single 
gust of predetermined magnitude and shape. Classical discrete models used by NASA 
during the Apollo and Space Shuttle Programs included a 9 m/sec quasi-square-wave gust 
with variable wavelength from 60 to 300 m, A later study derived discrete gust from a 
military specification (MIL-SPEC) document that used a “1-cosine” shape. The MIL- 
SPEC document contains a curve of non-dimensional gust magnitude as a function of 
non-dimensional gust half-wavelength based on the Dryden spectral model, but fails to 
list the equation necessary to reproduce the curve. Therefore, previous studies could only 
estimate a value of gust magnitude from the curve, or attempt to fit a function to it. This 
paper presents the development of the MIL-SPEC curve, and provides the necessary 
information to calculate discrete gust magnitudes as a function of both gust half- 
wavelength and the desired probability level of exceeding a specified gust magnitude. 
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1. INTRODUCTION 

Analysis of spacecraft vehicle responses to 
atmospheric wind gusts during flight is important in the 
establishment of vehicle design structural requirements 
and operational capability. Typically, wind gust models 
can be either a spectral type determined by a random 
process having a wide range of wavelengths, or a 
discrete type having a single gust of predetermined 
magnitude and shape. Classical discrete models used 
by NASA during the Apollo and Space Shuttle Programs 
included a 9 m/sec quasi-square-wave gust with 
variable wavelength from 60 to 300 m (NASA, 2000). A 
later study derived discrete gust from a military standard 
(MIL-STD) document that used a "l-cosine'' shape 
(Adelfang and Smith, 1998). The MIL-STD document 
contains a curve of non-dimensional gust magnitude as 
a function of non-dimensional gust half-wavelength 
based on the Dryden spectral model, but fails to list the 
equation necessary to reproduce the curve (DoD, 1990). 
Therefore, previous studies could only estimate a value 
of gust magnitude from the curve, or attempt to fit a 
function to it (Adelfang and Smith, 1998). Furthermore, 
the MIL-STD curve is based on a 1% risk gust 
magnitude 1 , so there was no way to determine gust 
magnitudes for, other risk levels. 

This paper presents the development of the MIL- 
STD curve, and provides the necessary information to 
calculate discrete gust magnitudes as a function of both 
gust half-wavelength and the desired probability levei of 
exceeding a, specified gust magnitude. The 
development background is an extension of the 
descriptions provided in Chalk et al, 1969, and Jones, 
1967. 

2. BACKGROUND 

For analyses of spacecraft vehicles, wind gusts can 
be treated as either random (spectral turbulence) or 
discrete. For random gusts, typical spectral models 
include the Von Karman and Dryden turbulence models. 
The Von Karman model has widely been considered the 
more “realistic” model when it comes to defining 
turbulence spectra. However, due to the computational 
complexity of the Von Karman model, the Dryden model 
is typically used in aerospace vehicle analyses. The 
longitudinal, lateral, and vertical Dryden spectra are: 
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where Q is spatial frequency (wavenumber), o is the 
turbulence standard deviation, L is the turbulence scale 
length, and the subscripts u, v, and w denote the 
longitudinal, lateral, and vertical components, 
respectively (NASA, 2000). The longitudinal component 
of turbulence is parallel to the steady-state wind vector, 
while the lateral and vertical components are 
perpendicular to it. The non-dimensional Dryden spectra 
are shown in Figure 1 . 

The spectra shown herein are in the spatial domain. 
The spectra can easily be transformed to the frequency 
domain by use of the Jacobian transform D = cu/V where 
« is radial frequency and V is the magnitude of the 
mean wind vector relative to the speed of the aerospace 
vehicle (NASA, 2000). 

An important detail that will be needed later is that 
the autocorrelation (lag correlation) is determined by 
taking the inverse Fourier transform of the Dryden 
spectra. The autocorrelation describes the correlation 
between gusts separated by a distance or time interval. 



Figure 1. Non-dimensional Dryden spectra for the 
longitudinal, lateral, and vertical components of 
turbulence. 




For the spatial domain, these autocorrelations are 
(Hogge, 2004): 
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where d is the lag distance. The autocorrelations as a 
function of non-dimensional values d/L are given in 
Figure 2. 
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Figure 3. Graphical depiction of the 1-cosine 
discrete gust. 


3. DEVELOPMENT OF DISCRETE GUST 
MAGNITUDES 


Generally, random gusts about a mean wind are 
considered to be normally distributed. This is true for 
each of the gust components (longitudinal, lateral, and 
vertical). Therefore, the random gusts, V, will have a 
probability density function (pdf) of the form: 
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where pv and ov are the mean and standard deviation of 
Figure 2. Autocorrelation as a function of non- the gusts. An initial gust, Vi, will be related to a gust 
dimensional values d/L. some distance away, V 2 , by the conditional probability: 


The discrete gust provides a “spike-type" input 
whose magnitude is based on information from the 
spectral model. Along with the gust magnitude, the gust 
gradient aver a specified temporal or spatial interval is 
also important. The classical shape of the discrete gust 
is that of a “1 -cosine" shape given in Equation 3 and 
shown in Figure 3 as a function of gust width (Chalk et 
al, 1969). 
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Here, V is the gust magnitude at distance x, and V m is 
the gust magnitude at d m , the gust half-width. Several 
values of can be chosen to tune the gust width to 
excite desired vehicle responses. The next section 
describes the development of the methodology to 
determine appropriate discrete gust magnitudes for 
given gust half-widths. 


/(F 2 |F,)=A^ (5) 

where f(Vi,Vg) is the joint pdf of the two gusts. The 
mean values of the gust, Vi and V?, are simply the value 
of the mean wind. The important factor in this 
development is the deviation of gusts around a mean, 
not the value of the mean itself. Therefore, these can be 
removed by setting them equal to zero. If the initial gust, 
Vi, is assumed to start at zero, then Equation 5 will 
reduce to: 
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where p is the correlation between Vi and V 2 . If we let 

e = (7) 




and substitute into Equation 6, we get the common form 
of the normal pdf (see Equation 4). 
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The cumulative distribution function (cdf) then becomes: 


deviation and length scale at this level are a = 7.72 m/s 
and L = 1230 m, The risk of exceeding the gust 
magnitude is 1%. Therefore, the probability to input into 
the inverse normal computer routine is 
P=1-(0.01/2)=0.995. The resulting gust magnitude is 
14.63 m/sec. Figure 4 shows the longitudinal gust 
magnitude for this example with gust half-widths in the 
range 1 to 10000 m. 
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Since Equation 9 can not be integrated in closed form, 
computer routines are relied upon to determine the 
normal cdf. To determine the gust magnitude value for 
V?, computer routines to calculate the inverse of the 
normal cdf are needed. Here, one would input the 
probability level (P), the mean of the gusts (jU - 0), and 
the standard deviation (e) To determine e, substitute 
the appropriate value of R from Equation 2 for p in 
Equation 7. 

Care must be taken when performing this 
caiculation. Recall that gust magnitudes are normally 
distributed and are two-sided. That is to say, the 
probability of having a gust magnitude greater than V 2 is 
the same as the probability of having a gust less than 
-Vz. Therefore, the probability level chosen must be the 
two-sided probability. Most computer routines compute 
the one-sided probability. To account for this, let 
P=1-{risfc/2), where risk is the desired probability of 
exceeding a given gust magnitude. 

Appropriate turbulence standard deviations, a, and 
length scales, L, are provided in Table 2-79b of NASA, 
2000 for various altitudes and turbulent cases (light, 
moderate, and severe}. Table 1 below lists the values 
from 1 to 10 kilometers for the severe case. 


Table 1. Severe turbulence standard deviations and 
length scales for the longitudinal <U), lateral (V), and 
ve rtical (W) turbulent components. 


Altitude 

(km) 

Standard Deviation 
(m/sec) 

Length Scale (m) 

U 

V,W 

U 

V, W 

1 

5,70 

4.67 

832 

624 

2 

5.80 

4.75 

902 

831 

4 

6.24 

5.13 

1040 

972 

6 

7.16 

5.69 

1040 

1010 

8 

7 59 

5.98 

1040 

980 

10 

7.72 

6.00 

1230 

1100 



Figure 4. Longitudinal gust magnitude for varying 
gust half-width, < 7 = 7.72 m/sec and L = 1230 m. 

A previous study (Adelfang and Smith, 1998) 
attempted to fit a function that reproduced the MIL-STD 
curve, which is non-dimensional gust half-width (dm/L) 
vs. non-dimensional gust magnitude (Wcr). The two 
curves are shown in Figure 5. Notice that the Adelfang 
and Smith curve is slightly conservative compared to the 
MIL-STD curve. Figure 6 is the same as Figure 4, but 
with the corresponding Adelfang and Smith curve 
included. The percent difference between the two 
curves is depicted in Figure 7. The Adelfang and Smith 
curve over estimate the MIL-STD curve by as much as 
13% in the 200 to 300 m gust half-width range. 

5. SUMMARY 

A development method of the MIL-STD non- 
dimensional discrete gust magnitude curve has be 
presented. The development allows for slight reduction 
in gust magnitudes determined by other studies 
(Adelfang and Smith, 1998). Also, the new technique 
allows for selection of various risk levels of exceeding 
gust magnitudes (see Figure 8). 
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As an example, a longitudinal gust magnitude for 
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Figure 5. Comparison of MIL-STD curve to Adelfang 
and Smith, 1998, curve. 



Figure 6. Same as Figure 4, but with Adelfang and 
Smith, 1998, curve included. 
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Figure 7. Percent difference between the Adelfang 
and Smith curve and the MIL-STD curve. 



Figure 8. Non-dimensional gust half-width (dm/L) vs. 
non-dimensional gust magnitude (Vm/o-) for various 
risk levels (MIL-STD Is the 1% curve). 
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Introduction 

Analysis of spacecraft vehicle responses to atmospheric wind gusts 
during flight is important in the establishment of vehicle design 
structural requirements and operational capability. Typically, wind 
gust models can be either a spectral type determined by a random 
process having a wide range of wavelengths, or a discrete type 
having a single gust, of predetermined magnitude and shape 
Classical discrete models used by NASA during the Apollo and Space 
Shuttle Programs included a o m/sec quasi-square-wave gust with 
variable wavelength from 60 to 300 m (NASA. 2000). A later study 
derived discrete gust from a military standard (MIL-STD) document 
that used a "i-cosine" shape (Adelfang and Smith, 1008). The MIL- 
STD document contains a curve of non-dimensional gust magnitude 
as a function of non-dimensional gust half-wavelength based on the 
Dryden spectral model, but fails to list the equation necessary to 
reproduce the curve (DoD. 1900) Therefore, previous studies could 
only estimate a value of gust magnitude from the curve, or attempt 
to fit a function to it (Adelfang and Smith. 1008). Furthermore, the 
MIL-STD curve is based on a i°e risk gust magnitude, so there was 
no way to determine gust magnitudes for other risk levels. The 
development of the MIL-STD curve is provided herein, and the 
necessary information to calculate discrete gust magnitudes as a 
function of both gust half-wavelength and the desired probability 
level of exceeding a specified gust magnitude are provided. 


Background 

For analyses of spacecraft vehicles, wind gusts can be treated as 
either random (spectral turbulence) or discrete For random gusts, 
typical spectral models include the Von Karman and Dryden 
turbulence models. The Von Karman model has widely been 
considered the more '‘realistic" model when it comes to defining 
turbulence spectra. However, due to the computational complexity 
of the Von Karman model, the Dryden model is typically used in 
aerospace vehicle analyses. The longitudinal, lateral, and vertical 
Dryden spectra are: 
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where S) is spatial frequency (wavenumber), a is the turbulence 
standard deviation, L is the turbulence scale length, and the 
subscripts u. v, and w denote the longitudinal, lateral, and vertical 
components, respectively (NASA. 2000) The longitudinal 
component of turbulence is parallel to the steady-state wind vector, 
while the lateral and vertical components are perpendicular to it. 
The non-dimensional Dryden spectra are shown in Figure 1. 



An important detail needed in the development is that the 
autocorrelation (lag correlation) is determined by taking the inverse 
Fourier transform of the Dryden spectra. The autocorrelation describes 
the correlation between gusts separated by a distance and are given by: 
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where d is the lag distance The autocorrelations as a function of non- 
dimensional values d L are given in Figure 2. 



Figure 3. Autocorrelation as a function of non- 
dimensional values d/L. 

The discrete gust provides a "spike-type - ' input whose magnitude is 
based on information from the spectral model. Along with the gust 
magnitude, the gust gradient over a specified temporal or spatial 
interval is also important The classical shape of the discrete gust is 
that of a "l-cosine" shape given in Equation 3 and shown in Figure 3 as 
a function of gust width (Chalk et al. 1060) 

V =0 x <0 

(3) 

v = 0 x>2 d m 

Here. V is the gusl magnitude at distance x, and V n is the gust 
magnitude at d„. the gust half-width. Several values of d_ can be 
chosen to tune the gust width to excite desired vehicle responses The 
next section describes the development of the methodology to 
determine appropriate discrete gust magnitudes for given gust half- 
widths 



Figure 3. Graphical depiction of the i-cosine discrete gust. 


Development of Discrete Gust Magnitudes 

Generally, random gusts about a mean wind are considered to be normally 
distributed This is true for each of the gust components (longitudinal, 
lateral, and vertical). Therefore, the random gusts. \\ will have a probability 
density function (pdf) of the form: 



where p v and o v are the mean and standard deviation of the gusts. An 
initial gust, 1',, will be related to a gust some distance away. by the 
conditional probability: 

( 5> 


where /0 W »* the joint pdf of the two gusts. The mean values of the gust 
V, and Vf. are simply the value of the mean wind. The important factor in 
this development is the deviation of gusts around a mean, not the value of 
the mean itself. Therefore, these can be removed by setting them equal to 
xero. If the initial gust. V v is assumed to start at zero, then Equation 5 will 
reduce to: 
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where p is the correlation between V, and V r If we let 

e = cr^l-p 2 (7) 

and substitute into Equation 6. we get the common form of the normal pdf 
(see Equation 4) : 
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The cumulative distribution function (cdf) then becomes: 
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Since Equation o can not be integrated in closed form, computer routines 
are relied upon to determine the normal cdf. To determine the gust 
magnitude value for Vg, computer routines to calculate the inverse of the 
normal cdf are needed. Here, one would input the probability level (P), the 
mean of the gusts (p = o). and the standard deviation (e). To determine e, 
substitute the appropriate value of R from Equation 2 for p in Equation 7. 
Appropriate turbulence standard deviations, o, and length scales. L. are 
provided in Table 2-7ob of NASA. 2000 for various altitudes and turbulent 
cases (light moderate, and severe). Table 1 below lists the values from 1 to 
to kilometers for the severe case. 


Table 1. Severe turbulence standard deviations and length scales for the 
longitudinal (l')> lateral (V), and vertical (W) turbulent components. 



Example Calculation 

As an example, a longitudinal gust magnitude for severe turbulence at the 
10 km level for a 500 m gust half-width is computed. The turbulence 
standard deviation and length scale at this level are a = 7.72 m/s and L = 
1230 m. The risk of exceeding the gust magnitude is iV Therefore, the 
probability to input into the inverse normal computer routine is 
P=i-(o.oi/2)=o.Q95 The resultinggust magnitude is 14.83 m sec. Figure 
4 shows the longitudinal gust magnitude for this example with gust half- 
widths in the range 1 to 10000 m. 
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Figure 4. Longitudinal gusl magnitude for varying 
gust half-width, a -7 72 m/secand t » 1230 ra 


Summary 

A development method of the MIL-STD non-dimensional discrete gust 
magnitude curve has be presented The development allows for slight 
reduction in gust magnitudes determined by other studies (Adelfang and 
Smith. 1008). Also, the new technique allows for selection of various risk 
levels of exceeding gust magnitudes (see Figure 5). 



Figure 4- Non-dimensional gust half-width (d„ L) vs. non-dimensional 
gusl magnltudeO'a, *) for various risk levels (MIL-STD Is Ihe iScurvc) 
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Figure 1. Non-dimensional Dryden spectra for the 
longitudinal, lateral, and vertical components of turbulence. 
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